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Abstract 

We complete the calculation of master integrals for massless three-loop form factors by 
computing the previously-unknown three diagrams with nine propagators in dimensional 
regularisation. Each of the integrals yields a six-fold Mellin-Barnes representation which we 
use to compute the coefficients of the Laurent expansion in e. Using Riemann £ functions of 
up to weight six, we give fully analytic results for one integral; for a second, analytic results 
for all but the finite term; for the third, analytic results for all but the last two coefficients in 
the Laurent expansion. The remaining coefficients are given numerically to sufficiently high 
accuracy for phenomenological applications. 



1 Introduction 



The quark form factor 7* — > and gluon form factor H ^ gg (effective coupling) are the 
simplest processes containing infrared divergences at higher orders in massless quantum field 
theory, and therefore are of particular interest in many aspects. They have, for instance, been 
used to predict the infrared pole structure of multi-leg amplitudes at a given order [1-4]. The 
form factors can also be exploited to extract resummation coefficients [5,6], and they enter the 
purely virtual corrections to a number of collider reactions (Drell-Yan process, Higgs production 
and decay, DIS). 

Besides phenomenological applications, a major motivation for obtaining analytic results at 
three-loop order and beyond is finding and understanding structures in massless gauge theories 
that generalize to an arbitrary number of loops. Much progress has been achieved recently in the 
prediction of all-order singularity structures in QCD [7,8], in conjectures about the all-orders 
behaviour of maximally supersymmetric Yang-Mills theories [9-17] and in investigations of the 
finiteness of N=8 supergravity [18-23]. It has also been shown, at two loops [24,25] and recently 
even for the matter contributions at three loops [26], that the soft anomalous dimension matrix 
in any massless gauge theory is proportional to the one-loop matrix. 

The two- loop corrections to the massless-quark [27] and gluon [28, 29] form factors were 
computed in dimensional regularisation with D = 4 — 2e to order e° and subsequently extended 
to all orders in e in ref. [30]. Two- loop corrections to order e° are also available for massive 
quarks [31]. The three-loop form factors to order e _1 (and e° for contributions involving fermion 
loops in the quark form factor) were computed in refs. [6,32]; see also ref. [33]. 

In order to calculate the quark and gluon form factors at higher orders in perturbation theory, 
the amplitudes are reduced to a small set of master integrals by means of algebraic reduction 
procedures [34-38]. At the three-loop level, the master integrals for massless form factors were 
identified in ref. [39], and results for certain subsets are available in the literature [34,39-42]. 
Among the three-loop master integrals, the genuine three-loop vertex functions are the most 
challenging ones from a computational point of view. They correspond to two-particle cuts of 
the master integrals for massless four-loop off-shell propagator integrals [43], which have been 
used in the calculation of the scalar i?-ratio [44]. In fact, such a correspondence, via two-particle 
cuts, between these two families of master integrals follows from the general result of ref. [45]. 
The derivation of the three-loop vertex integrals is of comparable complexity to massless four- 
loop propagator integrals. 

Working in dimensional regularisation and expanding the master integrals in a Laurent series 
in e, the finite part of the three-loop form factors requires the extraction of all coefficients through 
(polylogarithmic) weight sbH, i-e. coefficients containing terms up to 7r 6 or £|. 

Those genuine three-loop vertex functions which contain one-loop or two-loop propagator 
insertions were computed in ref. [39]. The three- loop master integrals which are sufficient in 
order to obtain the fermion loop contributions within a Feynman diagrammatic approach were 
computed in ref. [40] . The purpose of this Letter is to give the results of the remaining three dia- 
grams which have nine propagators each. We present analytic results for all but three coefficients 
in the e expansion, along with accurate numerical values for the remaining ones. 

*We prefer to use the term "weight" instead of "transcendentality" , because from a mathematical point of 
view there is no proof that £3 is a transcendental number. 
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Ag,l Ag t2 ^9,4 

Figure 1: Three-loop master integrals with nine massless propagators. The incoming momentum 
is q = pi + p2- Outgoing momenta are taken to be on-shell and massless, p\ = p\ = 0. 

2 Computational Methods and Results 

In this section, we list the results we obtained for the three- loop master integrals with nine 
propagators. They are depicted in Fig. [Hand labeled as in refs. [39,40]. All other diagrams with 
up to eight propagators were already given in the same references. We will work in dimensional 
regularisation with D = 4 — 2e. 

2.1 Diagram A 9 ;i 

The first diagram to be considered is ^4^1 which can be written as follows, 
d D k f d D l fd D r 

x \ (1) 

k 2 (k+ Pl ) 2 {k + lf (k-rf (l + r) 2 (l + p 2 ) 2 I 2 {r+ Pl f (r-p 2 f ' 

Here and in the following we tacitly assume that all propagators contain an infinitesimal +irj. 
The integral in eq. ([I]) can be written in terms of the following six-fold Mellin-Barnes (MB) 
representation [46-49]: 

dwB f dwp, 



„'j r 9 i -3-3e r°(l — e) f dwi f du>2 f dw 3 f dw4 f dw§ f 

A ^ = * 5r i- q ~ lT} \ tp^t J 2^ J 2^ J 2^ J 2^ J 2^ J 



C2 +i 00 

f dw 2 


Cg + i 00 

f dw 3 


c^ + i 00 

f dw^ 


J 2m 


J 2m 


J 2m 


C2 — i 00 


C3 — i 00 


C4 — i 00 



X 



r(-2e) J 2iri J 2m J 2m J 2m J 2m J 2m 

v(-w2) r(2 + w\ + w 2 ) r(-u) 3 ) r(ty 3 - w 2 - ^4) r(-u> 4 ) r(-^ 5 ) r(i + u>3 + ^5) 
r(— w 2 - 104) r(2 + ^3 + w 5 ) r(2 + w 5 + w e ) 

T(w 4: + w 5 - wi) T(-w e ) r(l + + u^) r(— 2 -3e - w 3 - w 5 ) 
X r(-l - 4e - w 5 ) T(-l - 3e - w 2 ) r(3 + e + wi + w 2 ) T(3 + 3e + w 5 ) 

x T(— 1 — 2e — W2 — W4) T(— 2 — 3e — UI5 — wq) T(2 + e + w\ + 102) r(u>i — u>5 — e) 
x r(3 + 2e + u> 2 + + ^5 + w 6 ) r(3 + 3e + w 3 + kj 5 ) T(3 + 3e + w 5 + u> 6 ) 
xr(-l- £ -wi)r(-l-e-w 2 ) , (2) 

where q 2 = (pi + p 2 ) 2 and 5 r = - — -. (3) 

(47r) ' r(l — e) 
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The representation ([5]) was obtained from eq. ([T]) by first introducing Feynman parameters; 
then, integrating over the momenta, loop by loop, and finally introducing MB parameters to 
decompose sums into products where appropriate. As usual, the contour integrals in the complex 
plane are along curves which separate left poles of T functions from right ones, where "left poles" 
are poles stemming from a T(. . . + w) dependence, while "right poles" stem from a T(. . . — w) 
dependence [49]. The most convenient choice for these contours are straight lines parallel to the 
imaginary axis, that is with constant real parts along the curves. According to refs. [47,48], 
these real parts, together with the parameter e, must be chosen in such a way as to have positive 
arguments in all occurring T functions in order to separate left and right poles in the desired 
way. In certain situations, such admissible straight contours and an appropriate starting value 
of e do not exist. This can be cured via the introduction of an auxiliary analytic regularisation 
(see, e.g., refs. [47,48]). 

The regularisation of the MB integral as described above, as well as the analytic continuation 
to e = 0, was done with the MB package [50] and, alternatively, with the MBresolve package [51], 
which is based on the strategy formulated in ref. [46]. (Within this latter strategy, straight 
lines for the contours along the imaginary axis are not required at the beginning.) These 
packages were also used for numerical cross checks. Moreover, we have also derived seven- 
fold MB representations in two different ways: using the AMBRE package [52] and starting from 
the MB representation of ref. [10] derived for general powers of the propagators for the tennis 
court diagram. The numerical evaluation based on these two MB representations was again 
performed with MB. In addition we performed numerical checks with the sector decomposition 
methods of [53, 54] and the FIESTA [55] package. (See ref. [56] for another implementation of 
sector decomposition.) 

As analytic techniques we apply Barnes's lemmas and the theorem of residues to the multiple 
Mellin-Barnes integrals, and insert integral representations of hyper geometric functions as well 
as ip functions and Euler's B function where appropriate. We also make use of the HPL [57-60], 
HypExp [61,62], and barnesroutines [63] packages, as well as an in-house implementation of 
the nested sums algorithm [64,65]. The final result for ^4g i reads, 



A 



9.1 



l Tj 



-3-3e 



1 | 1 

18^ + 2? + 



+ 



129 



■in* 
~3~ 



+ 



53 
18 
158 



4vr2 
~27~ 



+ 



<3 



/537 



578 



81 
322vr 4 



20tt 4 \ 1 



29 22vr 2 
T + 27 



6 

14 9a 238 , 

— T 2 Ca - -g-Cs 



2Cs 



+ 



2133 



405 3 

47T 2 + 158C 3 " — — " — TT 2 (3 + ^rC5 

135 3 3 



2398tt 6 
5103 



+0(e 2 



(4) 



We emphasize at this point that all terms in eq. ([!]) have been derived by purely analytic steps. 
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We also derived this result by evaluating first the following integral with a numerator: 



An) 
^9,1 



d D k f d D l fd 



(2n) D J (2n) D J (2ir) D 



k 2 (k+ Pl )' z {k + lf {k-rf (l + rf {1 + P2? I 2 {r + Pif (r - p 2 ) 



\2 ' 



(5) 



Starting from the above mentioned MB representation in ref. [10] for the tennis court diagram 
and setting one of the indices to minus one (the index corresponding to the numerator r 2 ), we 
obtained the following seven-dimensional MB representation: 



A 



(n) 
9.1 



ci +i oo C2+ioo c^ + ioo c^+i oo c^+ioo cg+ioo c<j-\-i oo 

-2-3e r 3 (l — e) f dw\ f dw2 f dw 3 f dwi f dw$ f dw$ f dw 7 



2iri J 2ttz J 2tti J 2iri J 2ni J 2iri 

c-^-i oo c 2~ 1 00 c^—ioo c^ — ioc c^—ioo eg— ioo cy — ioc 



r(-2e) J 2m 

c-^—ioc 

T(-W2)T(1 + w\ + w 2 + W3)T(-W3)T(-W4)T{1 + w\ + w 4 )r(-l - e — w% - w 3 ) 
r(l - w 2 )T{l - w 3 )T{l -2e + Wl + w 2 + w 3 ) 

T(-l — e — wi — w 2 — Wi) T(2 + e + w\ + w 2 + w 3 + w 4 ) T(-e + w\ + w 3 - w 5 ) 
r(-4e - w 5 ) r(l - w 6 ) T(l -w 4 - w 7 ) 

T(2 + 3e + w 5 ) r(— 1 — 3e — IU5 — u> 6 ) r(-io 6 ) T(l + io 5 + w 6 ) 

T(— 1 — 3e — Wi — w§ — W7) T(— e + wi + W2 — w$ — wq — w-j) T(—w§) T(—wj) 

r(— w\ + w 5 + 107) T(l + e - wi - w 2 - w 3 + w 5 + Wq + w 7 ) . 

i(«) 



(6) 



It turns out that at each order in the e expansion of Ag ( , the coefficients have homogeneous 
weight. This property turns out to be very helpful when one uses the so-called PSLQ algo- 
rithm [66] . Postulating that a given numerical result can be represented as a linear combination 
of certain constants (typically £ values and powers of ir) accompanied by rational coefficients, 
this algorithm solves for the latter. Starting from eq. ([6]) and using the MB and MBresolve pack- 
ages, we applied the Barnes lemmas whenever possible to MB integrals which appeared after 
expanding in e. At this point, at worst two-dimensional MB integrals were left. We calculated 
these integrals to an accuracy of 25 digits, which was sufficient to obtain very stable PSLQ 
results. For the term of highest weight we used the assumption that it is a linear combination, 
with rational coefficients, of 7r 6 and £f. This lead to the following result: 



An) 
^9,1 



-2-3 e 



1 



i rj 



7T- 



14C 3 47vr 4 



36e 6 
85 



18e 4 9e 3 



405e 2 



27 s 



1 



20C5 ) - + 



11607T 6 

5103 




+ 0(e) 



(7) 



Subsequently, we derived an analytic relation between Ag t \ and A 9 ( by means of a Laporta 
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reduction [35] using, independently, the AIR [37] and FIRE [38] packages. The result rea< 

(„) 8(2D - 7){2D - 5) (3D - 10) (3D - 8)(D - 3)(1311D 2 - 11764D + 26396) 
9,1 ~ 9(D - 4) 5 (3D - 14)(5D - 22) (g 2 ) 4 4 

80(2D - 7) (3D - 14) (3D - 8)(D - 3) 2 
+ 3(D - 4) 4 (5D - 22) (g 2 ) 3 5,1 

64(2D -7){D- 3) 3 (69D 2 - 580D + 1220) 

9(D - 4) 4 (3D - 14)(5D - 22) (g 2 ) 3 5,2 
8(3D - 14) (3D - 10) (3D - 8)(D - 3) 2 ^ (M ) 
(D-4) 4 (5D-22)(g 2 ) 3 5,1 
32(2D - 7)(D - 3) 3 (45D - 202) (M) 
3(D - 4) 4 (5D - 22) (g 2 ) 3 5,2 

64(2D - 7)(D - 3) 2 20(2D - 7)(5D - 18) 

+ 3(D - 4)(3D - 14)(5D - 22) (g 2 ) 2 6,1 9(D - 4) 2 (g 2 ) 2 6 ' 2 

8(2D - 7)(3D - 14)(3D - 10) (D -3) 
+ (D-4) 3 (5D-22)(g 2 ) 2 6,3 

2(3D - 14) _ (3D - 14) V 

(5D - 22)g 2 7 ' 3 2(D - 4) (5D - 22) 9,1 ' 1 ' 

The integrals A4, A^f\ and A^P ar e listed in the Appendix, while all other integrals with fewer 
than nine propagators are given in refs. [39,40]. For convenience, all the corresponding graphs 
are shown in the Appendix, Fig. [2j We have checked that eqs. dU),© satisfy this relation. 

2.2 Diagram Aj j2 

The next diagram we consider is Ag t 2- It reads 
. fd D k f d D l fd D r 

1 

X k 2 (k+ Pl ) 2 (k-l+ Pl ) 2 {k-r-lf {l + rf {I+P2? I 2 (r+pif (r - p 2 f ' 
Like Ag t \, it can be written in terms of a six-fold Mellin-Barnes representation. 

■ n3 r 2 ■ i-3-3,r 3 (i-e) 

i9,2 



Aoo = iSZ 



q 2 — it] 



3-3e r 3 (l — e) f dw\ f dw2 f dwz f dw^ f dwz, f dw, 



r(-2e) J 2-Ki J 2m J 2iri J 2m J 2iri J 2iri 

c-^—ioc c 2~ * 00 eg — ioo c 4~ * 00 c^— zoo eg — i 00 



r(— wi) r(2 + wi + u^) r(-u> 5 ) r(— w 2 + ^3 + ^4 + 1) r(w 5 - w 2 ) r(u> 5 - W4) 

T(l -w 4 + u> 5 ) T(-w\ -W3 + wq) T(l - w 2 + w 5 + w 6 ) T(2 - e + wi + w 3 + 104) 



5 We also thank Beat Todtli for correspondence on this point 



T(-w 6 ) T(w G + 1) rp 6 - u> 3 ) T(l + w 5 + w 6 ) T(-2 - 2e - wi - w 3 ) 

r(2 — 2e + w\ + u> 3 + 104) r(— 1 — u>i — 3e) T(l — e + W4 — IU5) T(3 + e + + w 2 ) 



x T(-l - e - w 2 ) T(l - e + w x + w 3 ) T(w 2 - w 4 - e) T(l + w 4 - e) T(-l 
x T(l - e + + w 3 + u> 4 - ui 5 - w 6 ) r(e - wi - u> 3 - u) 4 + u) 5 + w 6 ) 
x T(u;4 — W5 — e) r(2 + e + uii + u>2) T(3 + 2e + u;i + w; 3 + 104) . 

The techniques we apply are the same as before. The final result for Aq 2 reads 

-3-3e 



(10) 



.4. 



9.2 



iSf 



r 

r 2 

-9? 

+ ( 

+ 



- <7 — 1 rj 
5 



110 



+ 



9 3 
/170 16tt 2 



/ _ 20 

v T ~ "27 
4vr 2 166 



7tt 2 \ 1 



)? + ( 



50 _ 
373vr 4 x 1 

7^ 



177T 2 
^27~ 



91 



494 A 
+ ^C3 



1080 

1877T 4 



V 9 9 ' 9 540 
+ (-670.0785 ± 0.0326) + 0(e 



(11) 



The number for the finite term was obtained with MB [50]. Again, all pole terms in eq. (jlip have 
been derived by purely analytic steps. As in the previous case we performed an independent 
analytic calculation of an integral with a numerator, which again turns out to have homogeneous 
weight^ at each order in the e expansion: 

d D k f d D l fd D r 



A, 



(n) 
9.2 



(2vr) D J {2tt) u J {2tt) 



D 



(12) 



(I 'Pi? 



k 2 (k+ Pl ) 2 (k-l+ Pl ) 2 (k 



r — l) (l + r) (I+P2) I (r + pi) (r — p 2 ) 



2 ■ 



Again, we relate Ag j2 to Ag I by means of a Laporta reduction [37,38]. The result reads 



A 



(n) 
9,2 



16(2D - 7) (2D - 5) (3D - 10) (3D - 8)(D - 3)(70D - 303) 



3(D - 4) 4 (2D - 9)(5D - 22) (g 2 ) 4 



64(2D - 7) (3D - 8)(D - 3) 2 (19D - 84) 
3(D-4) 3 (2D-9)(5D-22) 
16(3D - 10)(3D - 8)(D - 3) 2 



7 2 ) 3 



A 



5,1 



64(2D-7)(D-3) ; 



A. 



5.2 



(D-4) 2 (2D-9)(5D-22) (g 2 ) 3 
8(2D - 7)(5D - 18)(16D - 71) 
3(D - 4) (2D - 9)(5D - 22) (q 2 ) 2 



AM) 
A 5,l + 



A 6 , 2 + 



3(D-4) 3 (2D-9) (g 2 ) 3 
64(2P-7)(3P-10)(P-3) 3 
3(D - 4) 3 (2D - 9)(5D - 22) (q 2 ) 3 
10(D - 3) 
(2D - 9)q 2 



A 



(M) 
5,2 



.4 



7.1 



8(D - 4) 2 



(2D - 9)(5D - 22)g 2 



^7,3 + 



8(3D - 13)(3P - 11) 
(2D-9)(5D-22)g 2 



[A 7A + A 7 , 



(3P - 14) 
(5P - 22) 



^8 - q 2 A 



9,2 



(13) 



'We thank Lance Dixon for the suggestion that this particular numerator generates homogeneous weights. 



In) 

We observed homogeneity of weights for ^4g 2 in all coefficients where we have an analytic result. 

To evaluate Ag n \ we started again with eq. (jlOp and repeatedly applied the Barnes lemmas 
as much as possible. The resulting integrals were treated numerically. Afterwards, the numbers 
were plugged into the right-hand side of eq. (I13p and an expansion in e was performed. Assuming 

(n) 

homogeneous weight helped to minimize the number of possible constants in A 9 2 • There is only 
one constant each at weight 0,2,3 and 4 (1, ir 2 , £3, 7r 4 ) and two constants each at weight 5 (tt 2 (3 
and Cs) arid 6 (Cf and 7r 6 ). Using PSLQ, we reproduced the result in eq. (fTTj) up to order 1/e 2 . 
We note that the MB integrals contributing to the finite part of Ag^ have dimensionality as high 
as five, and therefore prevent us from achieving an accuracy which is sufficient for a successful 
application of the PSLQ algorithm. 

For the integral with the numerator, the result reads 



A 



(n) 
9.2 



1 r\ 



-2-3 e 



2 7ir 2 91C 3 373vr 4 
9^ + 271 4 + "9^ + 1080? 



+ 



179 
~2T 



tt 2 C3 + 167C 5 ) - + (395.3405 ± 0.0326) + 0(e) 



(14) 



2.3 Diagram A94 

The last diagram we consider is Ag^. It reads 
;d D k f d D l fd D r 

^9,4 



(2tt) D J (2ir) D J (2tt) d 



(15) 



k 2 (k+pi) 2 (k — r) 2 (k — r — l) 2 (I + r) 2 (I + P2) 2 I 2 (r + p\) 2 (r— p 2 ) 2 
Like the previous integrals, it can be written in terms of a six-fold Mellin-Barnes representation, 



A 



9.1 



% 77 



-3-3e 



T 3 (l 



r(-2e)r(-l-4e) 



dw\ f du>2 f dii)3 f dw^ f dw^ f dwQ 



2iri J 2ni J 2iri J 2iri J 2ni J 2iri 

cj-ioo c 2~ * 00 eg — i 00 c^-z 00 c^ — ioc eg — i 00 



^ r(-u; 1 ) r(i + wi + u; 2 ) r(-w 3 ) r(i - w 1 + w 3 ) t(w 3 - w 2 ) r(i + w±) r(i + w 5 ) 

r(l — wi) T(w± + W2 — W3 — W4 + w$ — 2e) T(l — 2e + w\ + W2) 

r(— ^5) t(w4 -W5 + 1) t(w 5 - W4) r(-u> 6 ) r(i + u> 3 + ui 4 + w 6 - w 5 ) 
r(2 - to 1 + ^3 + w±) r(i - w 2 + ^3 + ^4 - ^5) r(2 + u>3 + u>4 + w e ) 

x r(— 2 — 3e — W4) T(w\ + — W3 — 2e) r(— wi — e) r(wi — e) T(— 1 — e — u>2) 

x r(— 2 — 3e — W3 — u>4 + — wq) r(l + e — wi — W2 + ^3 + ^4) 

x T(l + u>2 - e) r(2 + e + wi + u> 2 + w 6 ) T(3 + 3e + u> 3 + u> 4 + u>e) • (16) 
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The techniques we apply are the same as above. The final result for Ag^ reads, 

-3-3 e 



A 



9,1 



i r\ 



— ! — ( 



9e 6 9e 
+ fl7 + 



1 



10vr z 



27 



1 / 14 



47^ 
27 



12 U )1 



717T 2 

27 



200 



47vr 4 \ 1 



810 



+ (117.3999538 ± 0.0000032) - 

+ (1948.167043 ± 0.000025) + 0(e) 



(17) 



The two numbers were again obtained with MB [50]. All higher pole terms in eq. (|17p have again 
been derived by purely analytic steps. In the case of Ag^ a homogeneous- weight master integral 
also exists, with a numerator: 



^9,4 



d D k f d u l [d u r 



(2tt) d J (2?r)" J (2tt) 



(18) 



fe (ft+Pl) (k - r) (k-r-l) (I + r) (Z+P2) i (r+Pl) (?* — P2) 
For ^9,4, we also used the same alternative method of calculation as in the case of Ag t 2- This 
procedure was again based on a relation between A94 to Ag^, which follows from a Laporta 
reduction [37, 38] . The relation reads 

16(2D - 7) (2D - 5) (3D - 10) (3D - 8) (13 - 3)(348D 2 - 3037D + 6618) 



A 



(n) 
9,1 



9(D - 4) 4 (2D - 9) (3D - 14)(5D - 22) (g 2 ) 4 
32(2D - 7) (3D - 8)(D - 3) 2 (31D - 138) 



A, 



3(D - 4) 3 (2D - 9)(5D - 22) (q 2 ) 



2^ 



A, 



5,1 



128(2D - 7) (.D - 3) 3 (195D 2 - 1726D + 3816) 
+ 9(D - 4) 3 (2D - 9)(3D - 14)(5D - 22) (q 2 ) 3 5 ' 2 



128(2D - 7)(D - 3) s 



3(D - 4) (3D - 14)(5D - 22) (q 2 ) 2 ' 
8(2D - 7) (3D - 10) (D - 3)(7D - 30) 

(D-4) 2 (2D-9)(5D-22) (q 2 ) 2 
8(3D - 13)(3D - 11) 



+ 



8(2D - 7)(5D - 18)(28D - 123) 
9(D - 4) (2D - 9)(5D - 22) (q 2 ) 2 



A, 



(2D-9)(5D-22)g 2 



[A 7 ,4 + A 7j5 ) 



. 20(D-3) . 
A6 ' 3 + (2D - 9)q 2 A > 2 

(3D - 14) V 



2(2D - 9)(5D - 22) 



A 



9,4 



(19) 



For the integral with numerator, the result reads 



A, 



(n) 
9,4 



1 

9? 



' —17] 

Wtt 2 
27? " 



-2-3 e 



12C 3 47tt 4 



810e 2 



+ (206.7612077 ± 0.0000032) - 



(1237.300592 ± 0.000035) + 0(e) . (20) 



8 



As was the case for Ag 2 , the remaining Mellin-Barnes integrals are of too high dimension to 
allow for a stable PSLQ fit to determine the remaining coefficients analytically. However, we 
are confident that a dedicated effort to determine the remaining coefficients analytically will 
eventually be successful. 

3 Conclusions and Outlook 

In this Letter we have evaluated the three nine-propagator master integrals needed for computing 
the quark and gluon form factors to three-loop order. Each of the three integrals can be expressed 
in terms of a six-fold Mellin-Barnes representation from which we determine all coefficients 
through weight six in the Riemann £ function. One integral is given fully analytically up to 
order e. For the second one we give an analytic representation for all pole parts, and the third 
one is given analytically except for the coefficients of 1/e and e° in the Laurent expansion in e. 
The remaining Laurent coefficients are given numerically to an accuracy which is sufficient for 
all phenomenological applications. 

Note added: Our results for the coefficients of the three master integrals Ag 1, Ag^ and Ag^ 
partially overlap with those of ref. [67] where these integrals were evaluated in an indirect way. 
Agreement has been found for all common coefficients: terms up to e° for Ag^, up to 1/e for 
Ag t 2 and up to 1/e 2 for Ag^. To avoid confusion, we would like to point out that the convention 
for the overall prefactor used in ref. [67] differs from our one. For other recent progress on the 
three-loop quark form factor see ref. [68]. 
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A Additional integrals 

In this appendix we collect some additional integrals which appear in the Laporta reduction of 
the nine-propagator integrals with numerator. Three of the integrals are two-point functions of 
the MINCER type [42]. The fourth one is the eight-propagator butterfly graph. The diagrams 
are depicted in Fig. [21 together with diagrams which have been calculated in Ref. [40]. The 
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results for the two-point functions are 



Aa 



d D k r d u i r d u r 



{2ir) D J {2tt) 1j J {2ir) n (k + qf {k + lf (I + r) z r 



Di 



A 



(M) 
5,1 



iSl 


-Q 2 - 


d D k 


fd D l 



! 1] 



2-3 6 r 7 (l - e )r(-2 + 3e) 
r(4 - 4e) 



(21) 



(2tt) D J (2vr) w J (2ir) u k 2 (k + q) z {I + qf {I + rf r 



d D r 



\D ,2 



A 



(M) 
5,2 



d D k f d u l fd n r 



(2ir) D J (2tt) 1j J (2ir) n (r + qf k 2 I 2 (k + r) 2 (I + r) 



Di 



1-3 6 r 8 (l-e)r(e)r(-l + 2e) 
r(2 - 2e)r(3 - 3e) 



(22) 



2 7.2 ,2 



1-36 r 8 (l - e)r 2 (e)r(-l + 3e)r(2 - 3e) 



(23) 



T 2 (2 - 2e)r(2e)r(3 - 4e) 

The last integral which we give is the eight-propagator butterfly graph As ; b, depicted in Fig. [2 
It is not a master integral, but it is useful in quite a number of calculations. It is obtained from 
Ag t i by shrinking the horizontal propagator to a point. From the Laporta reduction we obtain 

■d D k f d D l fd D r 



(2tt) D J {2tt) u J {2tt) 



D 



k {k+piY {k-rf (l + rf (l+p 2 y I 2 (r+PiY (r - p 2 Y 



128(2D -7)(D -3) ; 



A5.2 + 



(D - 4) 3 (3D — 14) (q 2 ) 3 5 ' 2 ' 3(D-4) 2 (3D- U) (q 2 ) 2 
48(21) - 7) (2D - 5) (3D - 10) (3D - 8)(D - 3) , 

Aa . 



32(2D-7)(D-3) 5 



+ 



A>,i 



(D-4) 4 (3D-14) (q 2 ) A 
Hence the integral can be written entirely in terms of T functions. Its expansion reads 

n-2-3e 



(24) 



A 



8,B 



1 77 



1 

9e5 3e 4 



1 / 2vr 2 \ 1 / 



2n 2 



52 



: C3 



1 



2 . 52 M 7tt 4 n 1 



9 - 2tt 2 + — C3 - 

3 s 90 / e 

7vr 4 



+ (-27 + 67 r 2 -52 C 3 + ^ + fvr 2 C3-^C 



140 



+ (81 - 18vr 2 + 156C 3 
+0(e 2 ) 



111 
10 



68 



VCs + 140C 5 + 



2473vr e 
34020 



+ 



1136 



-a e 



(25) 
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